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\Q , Production and equilibration of quark-gluon plasma are studied within 

I the color flux-tube model, at the RHIC and LHC energies. Non-Abelian 

J> ' relativistic transport equations for quarks, antiquarks and gluons, are solved 

in the extended phase space which includes coordinates, momenta and color, 
(f-) . Before the chromoelectric field is formed, hard and semihard partons are 

produced via minijets which provide the initial conditions necessary to solve 
the transport equations. The model predicts that in spite of the vast difference 
between the RHIC and LHC incident energies, once the local equilibrium is 
reached, the energy densities, the number densities and the temperatures at 
the two machines may not be very different from each other. The minijet 
input significantly alters the evolution of the deconfined matter, unless the 
r~| ! color field is too strong. For the input parameters used here the equilibration 

time is estimated to be ~ 1 fm at RHIC and ~ 0.5 fm at LHC, measured 
from the instant when the two colliding nuclei have just passed through each 
other. The temperature at equilibration is found to be ~ 250 MeV at RHIC 
and ~ 300 MeV at LHC. 
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I. INTRODUCTION 



Lattice QCD studies indicate that the hadronic matter undergoes a phase transition 
to the deconfined quark-gluon plasma (QGP) phase, at a sufficiently high temperature (~ 
200 MeV) [Q. The energy density at which this transition is expected to occur is about 2 
GeV/fm 3 . Such a phase of deconfined partons did exist in the very early universe and is 
expected to be created in the laboratory, in ultrarelativistic heavy-ion collisions (URHIC). 
The Relativistic Heavy-Ion Collider (RHIC) at BNL (Au-Au collisions at (y/s) NN = 200 
GeV) and the Large Hadron Collider (LHC) at CERN (Pb-Pb collisions at (y/s) NN = 5.5 
TeV) will provide the best opportunity to study the formation of QGP. The main signatures 
of formation of QGP are sought in (1) J/ip suppression, (2) strangeness enhancement and 
(3) dilepton and photon productions ||. 

Different stages through which the evolution of deconfined quark-gluon matter proceeds 
are (1) pre-equilibrium (2) equilibrium and (3) hadronization. While the equilibrium stage 
is described by the well known hydrodynamic equations 0, it is a difficult problem to study 
the formation and pre-equilibrium evolution of QGP in URHIC. Nevertheless, the initial 
conditions which one assumes for the hydrodynamic expansion, can only be determined by 
a careful study of the pre-equilibrium stage. So it is necessary and important to study 
the pre-equilibrium stage, which would provide information on the equilibration time, and 
the bulk properties such as the temperature, the energy density and the number density of 
QGP. For example, the importance of the correct determination of the equilibration time, in 
connection with the J/ip suppression, is stressed by many people and it is a challenging 
task to determine this quantity accurately. 

In this paper, we study the pre-equilibrium evolution of QGP within the color flux-tube 
model. In this model, when two Lorentz-contracted nuclei collide head-on and pass through 
each other, they acquire a nonzero color charge (< Q >= 0, < Q 2 >^ 0), by a random 
exchange of soft gluons. The two receding nuclei act as color capacitor plates and produce 
a chromoelectric field between them [7]||. This chromoelectric field causes instability of the 
QCD vacuum and creates qq and gluon pairs via the Schwinger-like mechanism PHTTI. The 
partons so produced, collide with each other and also get accelerated due to the background 



field. These color charges also rotate in the color space according to Wong's equation [12 



We consider relativistic non-Abelian transport equations for quarks, antiquarks and glu- 
ons. These equations simultaneously include the convective flow terms, the background field 
acceleration terms, the qq and gluon production source terms, the collision terms and the 
terms which take into account rotations of these color charges due to the background chro- 
moelectric field. We solve these transport equations together with the Yang-Mills equations. 

As mentioned above, the initial chromoelectric field owes its existence to nonperturbative 
processes, viz., the soft-gluon exchanges that take place between the two nuclei during the 
collision. This model describes mostly the production of soft partons JI3UI3] . As pointed out 
in [15] and also in [14|, the (semi)hard processes (via perturbative QCD) take place before 
the soft-gluon exchanges. By adding these hard processes to the color flux-tube model, we 
attempt to cover the full range of particle production mechanisms, thereby including both 
soft and hard partons in the evolution of the deconfined parton matter. It has been shown 
by several authors that minijet production is dominant in the very early stages of the col- 
lision of two nuclei at the RHIC and LHC energies |L6]. Earlier studies [jlY , 13|, 14| of the 



2 



pre-equilibrium domain ignored the minijets. In this paper, we have calculated the minijet 
production at RHIC and LHC and have combined it with the color flux-tube model to study 
the production and evolution of QGP. A similar attempt was made by Eskola and Gyulassy 
15fl within the hydrodynamics framework which is applicable only after local equilibrium 



is attained. Here we study the pre-equilibrium evolution by solving relativistic non-Abelian 
transport equations with both the hard and soft partons taken into account. We address our- 
selves to the important question whether QGP, if formed, would attain substantially higher 
temperatures at LHC than at RHIC. We also investigate relative importance of the hard and 
soft production mechanisms. 

The plan of the paper is as follows. In section II, we describe our model with a detailed 
discussion of minijet production and of transport equations with their analytic solutions. 
Here we also present our computational procedure. In section III, we describe our numerical 
results and discuss them. We then compare our model with some other models of the early 
evolution of the matter in URHIC. Finally in section IV, we present our conclusions. 



II. MODEL 



In our model, all hard and semihard partons (i.e., minijets) are formed before the chro- 
moelectric field is created due to the exchange of soft gluons. In all earlier studies of pre- 
equilibrium evolution of QGP within the color flux-tube model |17 , 13| , jjfl , there was only 
the background field and no particles, at the initial time. In other words, the initial particle 
distribution functions were taken to be zero in all these calculations. In the present work, 
we obtain the initial particle distributions from minijets. With this initial condition, we 
then solve the transport equations at the RHIC and LHC energies. Our aim is to study the 
possible effect of the minijets on the equilibration process. 

The relativistic non-Abelian transport equations which we want to solve are presented 
in |Ji~8|Ji4}1 , in detail. In order to include the color charge in the phase space, one considers 
an extended one-particle phase space of dimension d = 6 + (N 2 — 1), (with N = 3). The 
extended phase space is taken to be the direct sum R & © G, where G is the (compact) space 
corresponding to the given gauge group. Thus, in addition to the usual 3 space coordinates 
and 3 momentum coordinates, one now has 8 coordinates corresponding to the 8 color 
charges in SU(3). In this extended phase space, a typical transport equation in the notation 
of reads as 



P^ + gQ a F"p»d£ + gf 



abc 



Q a Aydi 



f(x,p,Q) = C(x,p,Q) + S(x,p,Q). 



Here f(x,p, Q) is the single-particle distribution function in the extended phase space. The 
first term on the LHS of Eq. (fj) corresponds to the usual convective flow, the second term 
is the non-Abelian version of the Lorentz force term and the third term corresponds to the 
precession of the color charge, as described by Wong's equation [Q. C and S on the RHS 
of Eq. (]]]) are the collision and the source terms, respectively. Note that there are separate 
transport equations for quarks, antiquarks and gluons JT8[ . 



In order to solve the transport equations, we need the initial distribution functions fo : 
which we now obtain using minijet production cross sections. 
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A. Minijet production at RHIC and LHC 



Hard and semihard partons expected to be produced at RHIC and LHC are mostly 
minijets whose average momenta are not very large. One can calculate the minijet cross 
sections after fixing a momentum scale (pt) above which perturbative QCD (pQCD) is 
applicable. In the lowest order pQCD the inclusive (2 — > 2) minijet cross section per nucleon 
in AA collision is given by [|l5i|T6f 



a jet = J dp T dy 1 dy 2 '^?^Yl x i fi/A{xi,$) %i fj/A^Pr) crij^ki(s,i,u). (2) 

ijkl 



Here x\ and x 2 are the light-cone momentum fractions carried by the partons i and j from 
the projectile and the target, respectively, / are the bound-nucleon structure functions and 
?/i and y 2 are the rapidities of the scattered partons. The symbols with carets refer to the 
parton-parton cm. system. The o%j-^u is the elementary pQCD parton cross section. The 
required &ij-+u are given by 

4q, 2 _§ 2 _|_ d 2 

Pqq'^qq* = "^f [~ J 2 1 ( 3 ) 

Aa 2 i 2 + u 2 

VqH^'q' = -g-d^^J> ( 4 ) 



Aa 2 s 2 + u 2 s 2 +i 2 2 



2 



4a 2 s 2 + u 2 i 2 + u 2 2w 2 . 



^=M(P + tf)[^-^], (8) 



9a 2 ut us st, , . 

*»™ = #[3-|,- F -fid- (1°) 



Here a s is the strong coupling constant, q and g' denote distinct flavors of quark, and 
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s = x x x 2 s = Ap\ cosh 2 



2 yi -V2 



(ii) 



The rapidities yi, y 2 and the momentum fractions Xi, x 2 are related by, 

xi = p T (e Vl + e y2 )/^, x 2 = p T (e' yi + e~ y2 )/^. (12) 

The limits of integrations of rapidities yi and y 2 are given by, \y\ \ < ln( v /i/2pr+y 1 s/Ap 2 ^ — 1) 
and — hi(yfs/pT — e~ yi ) < y 2 < \ia(y^s / pr — e Vl ) . We multiply the above minijet cross 
sections by the well known factor K = 2, to account for the 0(a 3 ) contributions. 

The structure functions fij/A occurring in Eq. (|2]) are obtained from the model of Eskola 
et al. p|,|20|], and using the GRV HO 94 set of parton distributions for a free nucleon. 
To|, we choose minimum p? 



As argued in 
be applicable. 



2 GeV, above which pQCD is assumed to 



B. Initial conditions for plasma evolution from minijets 

For central collisions, the minijet cross section (Eq. (0)) can be related to the total 
number of partons (N) by 

N = T(0)a jet , (13) 

where T(0) = 9y4 2 /87ri? 2 4 is the nuclear geometrical factor for head-on AA collisions (for 
a nucleus with a sharp surface) ||16|| . Here Ra = 1.2 A 1 ^ 3 is the nuclear radius. The total 
transverse energy < > produced in nuclear collisions, from minijets, can be shown to 
be 

< E^ 1 >= T(0) f dp T p T dy x dy 2 ^-^^x x fi/A{.x u p 2 T ) x 2 fj/ A {.x 2 ,p 2 T ) a i:j _ k i(s,i,u). (14) 

J S ijkl 



We will use these relations (Eqs. (0), ( |T3"D and (0)) to determine the initial conditions which 
are required to solve the transport equations (Eq. ([]])) in the pre-equilibrium stage. 

In the color flux-tube model, we assume that the chromoelectric field is formed when the 
two nuclei completely pass through each other. The volume of the system at that instant 
is Vo = T^R\d, where d is the longitudinal thickness of the system. Note that d depends on 
the incident energy because of the Lorentz contraction and the presence of the wee partons 
in the incoming nuclei. We have taken d = 3 fm at RHIC for Au-Au collisions and d = 2 fm 
at LHC for Pb-Pb collisions El]. Using (Eqs. (g), (|T3|) and (0)), we calculate the initial 



number density (n ), initial energy density (eo) and initial distribution functions (/o) of 
quarks, antiquarks and gluons, which are needed to solve the transport equations. These 
are n = N/V , e =< E^ 1 > /V and f = dn /d 3 p. Here 

d 3 p = d 2 prdp z = pt d 2 pT cosh?/! dy x , (15) 

and /o can be easily extracted from Eq. ( |T3"D with the help of Eqs. (|2|) and flip]). 

In Table I we present our results for the initial conditions, obtained from minijet pro- 
duction at RHIC and LHC. 
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C. Solution of the relativistic transport equations 



In order to solve the transport equations (|I|), we need to specify the collision and the 
source terms occurring on the RHS. 

We measure the proper time r from the instant the two colliding nuclei have just passed 
through each other. In our model, hard partons (in the form of the minijets) are present in 
the system at r = and soft partons are created by the background field at later times. It 
is not possible to obtain the collision term from pQCD alone because the latter provides a 
framework to treat only the hard collisions; collisions among soft partons cannot be handled 
by pQCD. Including the effects of soft collisions, in the transport equations (1), from first 
principles (QCD) is a difficult problem. We make the relaxation-time approximation to 
simulate the effects of both soft and hard collisions and write the collision term as: 



C = -p^(f-D/r c 



(16) 



where w M is the four velocity, f eq is the equilibrium distribution function (Fermi-Dirac for 
quarks and antiquarks, and Bose-Einstein for gluons) and r c is the relaxation time. We take 
the same expressions for the source terms (S) as in [Tlfl . 

Following Bjorken's proposal ||, we express all physical observables in terms of the boost 
invariant quantities, namely 



(t 2 



2U/2 



) 1 , € = v-v> pt = (p -pi) 



-2\l/2 



(17) 



where 77 = tanh _1 (z/t) is the space-time rapidity and y = tanh _1 (^/p ) is the momentum 
rapidity. The transport equations (H) can be rewritten in terms of these variables. 
The formal solutions of the transport equations can be found to be, 



f q ,q,g( T > f > Pt, Q) = dr' exp 



■M.g\T , ,Z , ,PT,Q) f%g,g{T , ,£',PT,Q) 



Pt cosh £' 



+ exp(-r/r c )/ 9 



q,q,9 



where £' is given by 

f = sinh" 1 



T 



sinh; 



e + ^i f T dr" r" E{r") 
Ptt' Jt' 



(18) 



(19) 



and fo 9 ' 9 are the initial distribution functions of quarks, antiquarks and gluons obtained 
from minijets (as described in subsections IIA-B). The initial distribution function was 
assumed to be zero in all earlier studies of the color flux-tube model [17,13,14]. In the above 
equation, 9\ is the angle between the chromoelectric field and the color charge in the SU(3) 



color space ]14 



Since both the field and the particles are present simultaneously, we use the following 
relation for the conservation of the energy-momentum tensor: 



u H 1 mat "r u fi 1 f 



0. 



(20) 



Here 
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rpflU 

■*■ mat 



p^(3f q + 3f g + f g )drdQ, 



(21) 



and 



rpflU 



diag(E 2 /2; E 2 /2, E 2 /2, -E 2 /2). 



(22) 



In Eq. (pID, the factor 3 refers to the three flavors of quarks, dT = d 3 p/(2n) 3 pQ = 
PTdpTd£/(2ir) 2 , and dVt is the angular integral measure in the color space. We follow the 



same procedure as in ||14| , p^7| and finally obtain the evolution equation for the field as: 



roo roo r 

/ dp T p 2 T / sinhf / dn[3fq _ 3/ _ + fg] + ^ /Q )a\E(r)\^ 2 = 0. (23) 

JO J —oo J 

2 - 876 and C(5/2) = 1.342 is the 



rfr (2tt) 2 

Here a = a((5/2) exp(0.25/a), a = cig/2) 5 / 2 ^ 
Riemann zeta function. 

To solve Eq. (|23|) , we fix the form of the local temperature T(r), by demanding that the 
particle energy density e p (r) at any instant, differs negligibly from the equilibrium energy 
density at temperature T(r) at that instant. (It may be noted that T(r) occurs in Eq. (|18[) 
and hence in Eq. (1231), through / ei3 .) This allows us to write the temperature T(r) in terms 
of the particle energy density e p (r) as follows 



T(r) 



36e p ( r) 
5vr 6 



.1/4 



where 



rfr^(p%) 2 (3/ g + 3/ 9 - + / 9 ). 



(24) 



(25) 



We solve Eq. (p3| ) numerically to study time evolution of the chromoelectric field. 
Finally, the number density of the quark-gluon matter is given by 



n(r) = J dTdQp^u,(3f q + 3f q + f g ). 



(26) 



D. Computational procedure 

In the present work, we have determined E(t) and T(r) by imposing the following 



double self-consistency requirement. Starting with trial E{r) and T(r) and using Eqs. (18) 
we determine f q , f q and f g . These when substituted in Eq. (p5|) give us e p (r) which in turn 
gives a new temperature T(r) by Eq. (|24|). The new T(r) and the trial E(r) are again used 
in Eqs. fllSD , and with the new f q , f q and f g thus determined, the differential equation (|23|) 
is solved to get a new field E(r). This process is iterated until convergence is reached. This 
gives us a self-consistent set of E(r), T(r), e p (r), / 9 , f q and / 5 . 
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III. RESULTS AND DISCUSSION 



With the minijet initial conditions, we solve the transport equations (1) from the instant 
the two wounded nuclei have just crossed each other up to r ~ 1.5 fm. The solution of the 
transport equations allows us to determine the time evolution of the physical quantities such 
as energy density, number density and temperature. The determination of these quantities 
is important for the prediction of any signature of quark-gluon plasma. We now present our 
numerical results for these quantities. As in [ Hf , we take g = 4 and r c = 0.2 fm. For the 



initial field energy density e/(0), we take either or 300 GeV/fm 3 . This allows us to discuss 
the following three scenarios: (a) vanishing field, (b) e/(0) comparable to the initial minijet 
energy density eo, and (c) e/(0) ^> eo- Thus we can study the effect of the variation of the 
initial field on the calculated results. 

In all the figures, r = corresponds to the instant when the two colliding nuclei have 
just passed through each other, and e/(0) = 300 GeV/fm 3 unless stated otherwise. 

In Fig. 1, we have presented the parton energy densities for RHIC and LHC. It can be 
seen that the evolution of the parton energy density with minijet inputs at RHIC, is almost 
the same as that without the minijet inputs (dashed curve). This is because, at RHIC, the 
initial minijet energy density eo (— 3 GeV/fm 3 ) is negligible compared to e/(0); see Table I. 
At LHC, the two results are different because eo (— 178 GeV/fm 3 ) is comparable to e/(0). 
Thus the inclusion of the minijet input substantially alters the evolution of the parton energy 
density unless the chromoelectric field is too strong. 

We now discuss qualitatively the shapes of the curves in Fig. 1. It is easy to see from 
Eq. (|18|) that the parton energy density evolves as 

e p (r) =e(r)+e-^e . (27) 

Here, e(r) is the energy density if there is no minijet production, i.e., if the particle produc- 
tion is only from the field. The decrease in the energy density at earlier times (see the curve 
labelled LHC in Fig. 1) is due to the multiplicative factor e~ T ^ Tc occurring in Eq. (|27j). The 
further increase in the energy density is due to the fact that the particle production from the 
field starts dominating by this time (this is e(r) in Eq. (^)). We shall show below that the 
fall in the energy density after it reaches the maximum value is close to what is expected if 
the system has attained a local equilibrium and is expanding according to the hydrodynamic 
equations. Then it is clear from Fig. 1 that the equilibration of plasma is expected to occur 
earlier and the temperature at equilibration is expected to be somewhat higher at LHC. This 
will have a significant effect on all the predictions of QGP signatures. 

In Fig. 2, we compare the evolution of the parton energy density, at LHC, with and 
without the background field. This is necessary because the initial strength of the field 
cannot be calculated easily and as the incident energy increases the possibility of soft-gluon 
exchanges may become small. So we have presented our results for the two extreme cases (1) 
vanishing field energy density and (2) a very high initial field energy density (300 GeV/fm 3 ). 
The actual results are expected to fall within these two limits. Qualitatively similar results 
are expected at the RHIC energy too. 

Thus the background field representing the collective long-range effects plays an impor- 
tant role in the evolution of the system. If the minijet initial conditions based on pQCD 
are to play an important role at RHIC, then the field will have to be about two orders 



S 



of magnitude weaker than what we have assumed. In other words, if at RHIC, the initial 
field energy density e/(0) was not 300 GeV/fm 3 , but only 3 GeV/fm 3 , then e/(0) would 
be comparable to the initial minijet energy density eo(RHIC) ~ 3 GeV/fm 3 , and results 
qualitatively similar to those shown in Fig. 2 would be obtained at the RHIC energy too. 



In the HIJING ]22] and Parton Cascade models [21], the background field was completely 
neglected. 

Comparison of the two curves labelled 'Field' and 'Field+Minijets', in Fig (2), shows the 
effect of the minijets on e p (r). The memory of the minijet input is nearly wiped out at large 
t because of the exponential decay factor in Eq. (P7[). As a result, the large kinetic energy 
of the two colliding nuclei at LHC may not necessarily translate into large temperature of 
the equilibrated deconfined matter. Actual calculation shows that the temperature when the 
equilibrium is reached is ~ 300 MeV at LHC and ~ 250 MeV at RHIC, not very different 
from each other. Comparison of the two curves labelled 'Minijets' and 'Field+Minijets' 
shows the effect of the field on € p (t); this effect survives till relatively larger r. Note also 
that if there is no background field, the e p (r) decreases monotonically. 

The decay of the chromoelectric field is shown in Fig. 3. The decay is much faster 
at LHC than at RHIC, again suggesting an early equilibration at LHC. The evolution of 
the field strength has a significant impact on the acceleration of partons. If the chromo- 
electric field survives even after local equilibrium is reached, then the subsequent evolution 
of the system is governed by the chromoviscous hydrodynamic equations [fL5|] , instead of 
Bjorken's hydrodynamic equations ||. So it is very important to study the evolution of the 
chromoelectric field in the pre-equilibrium stage. We find that the equilibration time of the 
deconfined matter and the life-time of the field are nearly the same. 

In Fig. 4, we have presented the parton number densities of the system, for both the 
RHIC and LHC energies. The behavior of the results is the same as that of the parton 
energy densities in Fig. 1, because an equation similar to Eq. fl2"7|) holds for the number 
density too. The maximum number density at LHC is ~ 120 fm -3 , whereas at RHIC it 
is ~ 90 fm -3 . However, at large r , when the systems have equilibrated, the two number 
densities are nearly the same. These partons produce /affect all the signatures of QGP and 
hence their evolution plays a crucial role in QGP detection. 

We have fitted e p (r), n(r) and T(r), for r > 0.5 fm, at the LHC energy, with a func- 
tional form ar~ b , where a and b are constants. We find that b = 1.23, 0.70 and 0.31, 
respectively. According to Bjorken's hydrodynamics, these values are 1.33, 1 and 0.33, re- 
spectively. Comparing these two sets of values with each other, we conclude that the system 
has nearly equilibrated at r ~ 0.5 fm, at LHC. Similarly, at RHIC, the equilibration time 
is found to be ~ 1 fm. Recall that these time intervals are measured from the instant when 
the two colliding nuclei have just passed through each other. 

In Figs. 5 and 6, we have presented our results for the number densities of quarks 
plus antiquarks and gluons, respectively. Consider the curves labelled LHC in these two 
figures. The initial number density in Fig. 5 is much less than that in Fig. 6, because the 
minijet gluon production is much larger; see Eqs. (3)-(10). In Fig. 5, initially, the number 
density increases with time because of the large production of q and q from the field; see 
Eq. ( p6| ) which receives a separate contribution from each of the 3 flavors of quarks and 
antiquarks. Although the gluon production is 1.5 times more likely than that of quarks via 
the Schwinger-like mechanism [fTl]] , the counting in Eq. fl2"E| ) eventually builds up a larger 
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quark plus antiquark density than the gluon density. In Fig. 6, the LHC number density 
decreases with time because of the dominance of the exponential factor exp(— r/r c ) in Eq. 
([[ID, over the gluon production by the field. The energy densities for quarks plus antiquarks 
and gluons are plotted in Figs. 7 and 8, respectively. The behavior remains the same as in 
Figs. 5 and 6. 

We now compare our model with some well known models of the early stages of the ultra- 
relativistic heavy-ion collisions, namely the parton cascade model (PCM) [T2T|], the HIJING 
model p2| , the McLerran-Venugopalan model (MVM) p3fl , etc. 

As is well known, PCM is a pQCD-based model. It neglects the long-range collective 
effects which we have attempted to incorporate here by means of a background field. The 
PCM starts by defining a classical phase-space distribution function for the two incom- 
ing nuclei, which is then evolved by means of the Boltzmann equation retaining only the 
convective-flow term and the collision term in Eq. (1). Now it is known that there are wee 
partons in each of the two nuclei and each wee parton spreads over the entire width of the 
nucleus |p4 |. It is difficult to define a classical distribution function for such a system of 
particles. Indeed the incoming nuclei are in pure quantum mechanical states. Hence solving 
the classical transport equation right from the instant the two incoming nuclei start touching 
each other, is a questionable procedure. 

A classical distribution function can be more justifiably defined when the spatial extent 
and the time intervals are larger than the average de Broglie wave length of the partons. 
In our model, we have started the classical evolution at a time when the two nuclei have 
completely crossed each other. At that instant, in our model, there are only hard partons 
(Pt > 2 GeV/c); the soft partons are created later when the system has grown even bigger. 
The initial condition for the hard part is obtained via the minijet calculation where no 
approximation is involved. The initial condition for the soft part is taken through the 
creation of a coherent chromoelectric field. 

The HIJING model too is a pQCD-based model. It combines pQCD processes with string 
phenomenology for non-perturbative soft processes. It assumes that the parton distributions 
can be approximated by thermal phase-space distributions with non-equilibrium fugacities. 
Thus the issue of thermalization of the deconfined matter in the framework of the kinetic 
theory is not addressed. 

Thus we have presented a model which is different from PCM and the HIJING Model. 
We believe that this model has some desirable features. For simplicity, we have neglected the 
interactions among the minijets before r = 0. Our main aim has been to see how the minijet 
input modifies our earlier analysis JT4j. It will be interesting to consider these interactions. 
Although pQCD provides a general framework for this purpose, including these interactions 
is quite nontrivial. 

Our approach has some similarities with the McLerran-Venugopalan model. The MVM 
is a classical effective field theory description of the small- x modes in very large nuclei at very 
high energies. This effective theory contains a scale fi which is proportional to the large gluon 
density at small x. The large gluon density ensures that even if the coupling is weak, the fields 
may be highly nonperturbative. They argued that the classical fields corresponding to the 
saddle point solutions of the effective theory are the non-Abelian analogue of the Weizsacker- 
Williams fields in classical electrodynamics. Now in a collision of such Weizsacker- Williams 
fields, the production of incoherent partons with transverse momenta p T » a s jj, can be 
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handled by means of pQCD and their propagation in the coherent field can be studied with 
the help of a classical transport equation [25]. In our model too both incoherent hard partons 
and coherent field (soft part) are taken into account. So in this sense our approach is similar 
to MVM. 

Real-time simulation of the full, nonperturbative evolution of the classical non-Abelian 
Weizsacker- Williams fields on lattice, for the gauge group SU(2), has recently been studied 
by Krasnitz and Venugopalan PB|; see also |27|. The calculation in [EHJ incorporates coher- 



ence effects which become important at small x and small p?, while reproducing simultane- 
ously the standard minijet results at large p?- However, the important issue of equilibration 
of gluons was not addressed. It will also be interesting to see both the coherent field and 
incoherent partons simultaneously taken into account in such a simulation. 



IV. CONCLUSIONS AND FUTURE PROSPECTS 

We have studied the production and equilibration of the quark-gluon plasma expected to 
be formed in ultrarelativistic heavy-ion collisions, at RHIC and LHC, within the color flux- 
tube model. The distribution functions of quarks, antiquarks and gluons are defined in the 
extended phase space of dimension 14 in SU(3). We have solved the non-Abelian relativistic 
transport equations for these distribution functions. The quarks, antiquarks and gluons are 
produced from the background field which is formed due to soft-gluon exchanges, via the 
Schwinger-like mechanism. The initial distribution functions of partons are calculated from 
the minijet production at RHIC and LHC. 

The model predicts that 

(i) The inclusion of the minijet input substantially alters the evolution of the bulk properties 
of the deconfined matter unless the chromoelectric field is too strong. 

(ii) Equilibration of the deconfined matter is expected to be faster at LHC than at RHIC. 
Estimated equilibration times are given. 

(ra) Surprisingly, in spite of the large difference between the LHC and RHIC colliding 
beam energies, the differences between the temperatures and equilibrium energy and number 
densities attained at the two machines, may not be very large. This is because the memory 
of the minijet input decays exponentially with time. Estimated numerical values for these 
physical quantities are given. 

(iv) Unlike the minijets, the effect of the background field survives over a longer time. 

(v) In the absence of a background field, the parton energy density and the number density 
decrease monotonically with time. 

(vi) The equilibration time of the deconfined matter and the life-time of the chromoelectric 
field are nearly the same. 

In summary, we have combined the non-perturbative aspects of QGP (background chro- 
moelectric field formation) with the perturbative aspects (minijet production) and have 
studied the evolution of the quark-gluon plasma. 

We have used in this work the Schwinger-like source terms p|-4TT| for parton produc- 
tion from a chromoelectric field. This mechanism is strictly applicable only for a constant, 
uniform field. However, in a heavy-ion collision, as the system evolves, the field acquires 
space-time dependence. Particle production in such a field has been studied and an appro- 
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priate source term has been derived in |28] , which will be incorporated in the pre-equilibrium 
evolution of the QGP and results will be reported elsewhere. 



Note added: 

Recently Bloch et al. have studied the thermodynamics of strong-field plasmas [[HJ. In 
the present paper, we used a classical transport equation while they have used a quantum 
equation. Our equation is non-Abelian, theirs is Abelian. Our initial conditions are different 
from theirs. Our source term is also different from theirs. But still the time dependence of 
the various quantities (their Figs. 1-3) is qualitatively similar to that obtained by us. This 
shows the robustness of our findings. 



ACKNOWLEDGMENTS 



We thank R.V. Gavai for many useful discussions. 



12 



REFERENCES 



[1] A. Polyakov, Phys. Lett. B 72, 477 (1978); L. Susskind, Phys. Rev. D 20, 2610 (1979); 
C. Borgs and E. Seiler, Nucl. Phys. B 215 [FS7], 125 (1983); Comm. Math. Phys. 91, 
329 (1983). For the status of the lattice QCD, see Proc. of the 17th Intl. Sympo. on 
Lattice Field Theory (Lattice'99), Pisa, Italy, June- July 1999, to appear in Nucl. Phys. 
B (Proc. Suppl.). 

[2] For the status of the various signatures, see the Proc. of the 14th Intl. Conf. on Ultra- 
relativistic Nucleus-Nucleus Collisions (QM'99), Torino, Italy, May 1999, to appear in 
Nucl. Phys. A. 

[3] J.D. Bjorken, Phys. Rev. D 27, 140 (1983). 

[4] G.C. Nayak, JHEP 02 (1998) 005. 

[5] G.C. Nayak, Phys. Lett. B 442, 427 (1998). 

[6] X.-M. Xu, D. Kharzeev, H. Satz and X.-N. Wang, Phys. Rev. C 53, 3051 (1996). 

[7] F.E. Low, Phys. Rev. D 12, 163 (1975). 

[8] S. Nussinov, Phys. Rev. Lett. 34, 1286 (1975). 

[9] J. Schwinger, Phys. Rev. 82, 664 (1951). 
[10] A. Casher, H. Neuberger and S. Nussinov, Phys. Rev. D 20, 179 (1979). 
[11] M. Gyulassy and A. Iwazaki, Phys. Lett. B 165, 157 (1985). 
[12] S.K. Wong, Nuovo Cimento A 65, 689 (1970). 
[13] G.C. Nayak and V. Ravishankar, Phys. Rev. D 55, 6877 (1997). 
[14] G.C. Nayak and V. Ravishankar, Phys. Rev. C 58, 356 (1998). 
[15] K.J. Eskola and M. Gyulassy, Phys. Rev. C 47, 2329 (1993). 

[16] K.J. Eskola, K. Kajantie and J. Lindfors, Nucl. Phys. B 323, 37 (1989); K.J. Eskola 
and K. Kajantie, Z. Phys. C 75, 515 (1997); X.-N. Wang, Phys. Rev. D 43, 104 (1991). 

[17] K. Kajantie and T. Matsui, Phys. Lett. B 164, 373 (1985); B. Banerjee, R.S. Bhalerao 
and V. Ravishankar, Phys. Lett. B 224, 16 (1989); A. Bialas, W. Czyz, A. Dyrek and 
W. Florkowski, Nucl. Phys. B 296, 611 (1988). 

[18] H.-T. Elze and U. Heinz, Phys. Rep. 183, 81 (1989). 

[19] K.J. Eskola, V.J. Kolhinen and PV. Ruuskanen, Nucl. Phys. B 535, 351 (1998). 
[20] K.J. Eskola, V.J. Kolhinen and C.A. Salgado, e-print |hep-ph/9807297 . 
[21] K. Geiger, Phys. Rep. 258, 237 (1995). 
[22] X.-N. Wang, Phys. Rep. 280, 287 (1997). 

[23] L. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 (1994); D 49, 3352 (1994). 
[24] J.D. Bjorken in Current Induced Reactions: Proc. of the 7th Intl. Summer Inst, on 

Theoretical Particle Physics, Hamburg 1975, edited by J.G. Korner, G. Kramer and D. 

Schildknecht, Springer- Verlag, Berlin, 1976, (Lecture Notes in Physics, Vol. 56), page 

93; J.D. Bjorken in Proc. of the Summer Inst, on Particle Physics, Stanford, July 1973, 

edited by D.W.G.S. Leith and S.D. Drell, Vol. 1, page 1. 
[25] L. McLerran, private communication. 

[26] A. Krasnitz and R. Venugopalan, Nucl. Phys. B 557, 237 (1999). 

[27] W. Poeschl and B. Mueller Phys. Rev. D 60, 114505 (1999); S.A. Bass, B. Mueller and 

W. Poeschl, J. Phys. G 25, L109 (1999); A. Kovner, L. McLerran and H. Weigert, Phys. 

Rev. D 52, 3809 (1995); D 52, 6231 (1995). 
[28] R.S. Bhalerao and V. Ravishankar, Phys. Lett. B 409, 38 (1997). 



[29] J.C.R. Bloch, CD. Roberts and S.M. Schmidt, e-print frucl-th/ 99 10073 



13 



TABLE I. Initial conditions for pre-equilibrium evolution of QGP, obtained from the minijet 
calculation. 





(^) NN (GeV) 


N 


V (fm 3 ) 


n (fm 3 ) 


e (GeV/fm 3 ) 


RHIC 


200 


504 


459 


1.1 


3 


LHC 


5500 


17422 


318 


55 


178 



Figure captions 

FIG. 1. Parton energy density vs proper time. The dashed line refers to the calculation 
without the minijet s. 

FIG. 2. Parton energy density vs proper time, at the LHC energy. The curve labelled 
Field (Minijets) is based on the calculation without the minijets (background field). The 
curve labelled Field + Minijets is based on the full calculation with both the minijets and 
the background field taken into account. 

FIG. 3. Chromoelectric field vs proper time. The dashed line as in Fig. 1. 

FIG. 4. Number density vs proper time. The dashed line as in Fig. 1. 

FIG. 5. Number density for quarks and antiquarks vs proper time. The dashed line as 
in Fig. 1. 

FIG. 6. Number density for gluons vs proper time. The dashed line as in Fig. 1. 
FIG. 7. Energy density for quarks and antiquarks vs proper time. The dashed line as in 
Fig. 1. 

FIG. 8. Energy density for gluons vs proper time. The dashed line as in Fig. 1. 
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